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ABSTRACT 



This thesis is the result of literary research on the subject of 
sensitivity in linear feedback control systems. It is a synopsis of 
information obtained from various technical publications and is designed 
to give the reader information on the theory of sensitivity and its 
application in design problems. Included are 

(1) Definition of sensitivity function and root sensitivity 
for incremental variations of single parameter and their 
design application. 

(2) Definition of sensitivity function applied to the large 
parameter variations and its design application. 

(3) Special topics such as specification for sensitivity 
function, sensitivity integrals, and introduction to 

<* ' < i 

multi-parameter sensitivity. 
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I 

INTRODUCTION 

One of the foremost properties of feedback is its ability to 
reduce the sensitivity of a system to variations of the system para- 
meters. Suppose there is an element in the system which must be used 
because of its special properties. However the element is sensitive 
to environmental conditions (such as temperature and pressure), or is 
subjected to aging, or has wide manufacturing tolerances such that when 
it is replaced the parameters of the new element may be markedly diff- 
erent from those of the old element. Whatever the reason may be, 
suppose the values of the element parameters vary significantly. It 
is desired that despite these variations certain system properties, 
such as the input-output transfer function or output impedance should 
remain substantially constant. The troublesome element could perhaps 
be replaced by a better one, but this may be very expensive or imposs- 
ible. In any case we would like at least to have an alternative to 
replacing the element. In most cases the troublesome elements are the 
active elements, vacuum tubes, transistors, energy conversion devices 
such as motors and generators. But they may also be passive elements 
or transducers. In general the system assumed in an analytical study 
will never precisely match the actual physical system: it is therefore 
important to know the effects of possible variations. 

A general knowledge of the effects of parameter changes can be 
used as a guide to system modifications which would improve the over- 
all performance. For all these reasons measures of closed loop system 
sensitivity to open loop system parameter variations are an integral 
part of the analysis problem. 
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Notions about system sensitivity were in the forefront when the 
feedback concept was Initially developed. This was natural, even 
unavoidable, since feedback systems possess the fundamental physical 
property that the effects of variation in the forward loop, whether they 
are taken as changes in open loop transfer function G(s), or as depart- 
ures from strict linearity or from freedom from extraneous noise, are 
reduced by the factor 1 + G in comparison with the effects which would 
be observed in a non feedback system. Accordingly, sensitivity measures 
were indispensable to any rational discussion of feedback systems and 
a useful, classical definition of sensitivity was made. Except for a 
minor modification the definition of sensitivity, given by Bode, remained 
unchanged for over a decade. Perhaps this static nature surrounded by 
dynamic growth in most other areas of feedback systems engineering, 
made the concept fade with little exposure to sensitivity concepts. 

In recent years classic sensitivity has become a more popular subject 
in automatic control and also in network synthesis. Finally, the 
emphasis on pole-zero specifications for system characteristics gave 
rise to new concepts of sensitivity, with associated new measures. 

These measures called here "gain," "open-loop pole," and "open loop 
zero" sensitivities, were evolved to account for changes in the posi- 
tion of closed loop poles due to shifts or changes in open loop gain, 
poles or zero. 

This thesis is the result of a literary research on the subject 
of sensitivity in linear feedback control systems. In section II, the 
system sensitivity and root sensitivity, for a system subjected to 
incremental parameter variation, are discussed. Based on the meanlns 
of open loop pole and zero sensitivity, a design technique for 
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compensating the system is also introduced. In Section III, the large 
parameter variation and a new definition of sensitivity is presented. 

It is shown that when the leakage transmission between input and output 
is negligible, the new definition and the classical definition of the 
system sensitivity are identical. In Section III, an attempt is made 
to utilize Nichols Charts for evaluation of the sensitivity at a part- 
icular frequency range. A design philosophy and technique is also dis- 
cussed. Section IV discusses some special topics which are of parti- 
cular interest today and could be selected as interesting topics for 
further investigation. 
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SENSITIVITY AND INCREMENTAL PARAMETER CHANGES 
II-I - System Sensitivity and Root Sensitivity 

Several definitions of sensitivity have been offered in the past. 
Bode defined the sensitivity of the overall transfer function T to the 
plant parameter k as 

-T _ }J»k k 

k bT 

T 

Horowitz used the inverse of Bode^s definition, 

bT 

-T \ jUT T 

k " i Av k bk 

k 

T 

Defined one way or another, S^ is generally known as "classical sensi- 
tivity" or "system sensitivity," since it involves the system transfer 
function T. Several other definitions have also been suggested, and 
each one of these definitions concentrates on some particular engineer- 
ing or mathematical simplicity. A. Schulke suggests since sensitivity 
relates changes in the transfer function of a system with respect to 
changes in parameter of the system, then one can define, 

1 _ _1_ _ b T 

sensitivity S " ^ /n,x 

and in general T is a complex and may be expressed as 

T(x) - A(x) + jB(x) 

where A(x) is the system attenuation function and B(x) is the system 
phase function, and from that, 

4" ' »i!*r ■ ~4r ~ f A(x > + J B <*>} 

X ^ 

- « 
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Then assuming that x is real , it is apparent that the real part of 
| is directly related to the system attenuation and that the imaginary 
part, is directly related to the system phase function. The use of 
S * which Schulke calls it "sensitiveness" permits comparison of both 
attenuation and phase functions relative to changes in a common para- 
meter. 

The most adopted definition of system sensitivity these days is 
the one suggested by Horowitz, i. e. 



g ^ 3 hJa. u, ^ 

k IU. k 



and we adopt it through the coming discussions. 

In recent years with the increased utilization of the pole-zero 
approach, it has become increasingly important to examine variations 
in position of the poles and zeros of the network function due to 
changes in network parameters. This is known as the "root sensitivity." 
Formal definitions of root sensitivity vary from author to author. 

Some define it as, 



Jkx 



b X 



where (~ s j) * s the root of the characteristic equation and x a para- 
meter which could be the gain K, or a pole or a zero . Huang uses the 



definition 



b s. 



s 



X 



X 
8 . 



S x 



Me Euer aid Stapleford prefer different definitions for sensitivity to 

S 1 S< l F J 

gain, , and se^sitivitv to poles, Sp' . or sensitivity to zeros, . 



J 



Again each one of these definitions leads to a particular engineering 
or mathematical simplicity utilized by that particular author. 

II-2 - Sensitivity Relation in Control Feedback Systems 

It has been noted by Bode that for any linear system the transfer 
function T(s,K), relating a response to an input function can be in 
the bilinear form of, 



and 



T(s,K) - 



S T _ M fcl 
K 



A(s) + KB(s) 
C(s) + KD(s) 

u 

- T 

JlL. 

K 



Both numerator and denominator of T are of the form, 

H(s,K) - q(s) + Kp(s) 

where K is a variable parameter and p(s) and q(s) are polynomials. The 
roots of H(s,K) » 0, which characterize H, depend on the value of K. 

H(s,K) - 0 constitutes an implicit relationship between s and K, and 

as K takes on all positive values, a number of curves in the S-plane 
are obtained. These curves are the loci of the roots of H called the 
"root loci." 

The root locus was introduced by Evans in 1948. He used it main- 
ly for the investigation of control systems with a transfer function 
T, in the form of, 

T . G Kp(b) 

1 + G q(s) + Kp(s) 

where G « Kp/q is called the forward gain in a unity feedback con- 
trol system. For such systems, only the denominator is of interest. 

As the result, there arose the generally accepted terminology in which 
the zeros of p and q were called open loop "poles" and "zeros," 
respectively and those of 1 + G ■ 0 are called "closed loop poles." 
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II-3 - Relationship Between System Sensitivity and Pole Zero 
Sensitivities 



The transfer function T(s) can generally be written in the form 



of 



T<s) 



S l rr (s +Z) 
m 

TT(g + 8 ) 
1 3 



(2-1) 



and let Z y s ^ and g be functions of some parameter x (such as temper- 



ature, gain of some element, etc.). Then 

n m 

U’l » + Sms + z ) - 2 Ms + s ) 

i 3 i 3 



and differentiating above we get 



dT 



dx 



n \ 7 

i & + 2 III 1 

g « 1 



m 






'ix s + Z. 



U- JL 

1 dx s + 



dividing both sides by dx 

x 



S 



dT 

T 

dx 

x 






n 






rf- + S * n 

^ 1 dx s + Z. 



m s 

- Zi * _1 



dx S + 8. 



(2-2) 



x J J 

Equation (2-2) gives the relationship between the classical definition 
of sensitivity and the closed loop pole and zero sensitivity, where 
the closed loop pole and zero sensitivities are defined as 



ds 

* -L 



❖ 



dx 

d Z 
dx 



(2-3) 



(2-4) 



This definition is not necessarily the most useful definition at all 
situations. 

Pole Sensitivity (closed loop pole): 

The root locus is defined by: 
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H(*,K) - q(s) + Kp(s) 



(2-5) 



hence- 



k - _ q(s) 

p(s) 



TT (s + p ) 

. L 

TT(S + Z ) 
i-1 



(2-6) 



and 



n m 

JU [ R ] = 2Ai(s + p.) - 2/^(8 + Z.) +Tj 
1=1 J 1=1 



(2-7) 



and 



& » da 
K 



n i tn i 

j5l + Pj^ i5 + Z i) 



Also since 



(2-8) 



K - - _g.(s> 

P(s) 



then 



. . q*p - p'q 

ds 2 

P 



(2-9) 



and 



dK 

ds 






£l ‘ El 

q p 



(2-10) 



For simplicity of notation Equations (2-6) and (2-8) can be written as. 



K - - TT (8 + s t ) i 



dK . 

- - ds 



i s + s. 



(2-11) 



( 2 - 12 ) 



which has to be evaluated at s = -s ^ . The roots of H ■ 0 are 

denoted as ( _s j) to distinguish them from the variable s which may of 

course assume any value. In other words -s is a point on the root locus. 

J 



16 



And 



if - is a root of q(s). 



if - is a root of p(s). 



i - - 1 

From Equations (2-10) and (2-12) we can write, 
8 . 



j 



ds 



dK 2 

K i ■ 



S + 8, 



IE1 



qp' - pq' 



S"8 , 



IJL 



q' + Kp' 



(2-13) 



where is the closed loop pole sensitivity with respect to K and 

sometimes referred to as just pole sensitivity with respect to K. 

Furthermore the residue of a function A(s) at pole s*a is defined as 

B(s) 



Residue 
at s= a 



J&il. 

B(s) 



lim (s - a) A(s) 
s-*a E(s) 



« lim (s - a ) A* (s) + A(s) 
s-*a B'(s) 



lim A(s) 
s-»a B'(s) 



which is true for simple poles. Noting that 

^s 

applying above residue theorem we get 



q' + Kp' and 



j 



’j 



K 



_KP_ 



q' + Kp' 



where Q is the residue of 



KJL 



> 8 J 



-Q, 



(2-14) 



T(s) 



K P( S > 



q(s) + Kp(s) 



at point s ■ -Sj. The same result can be obtained directly from (2-5) 



4 s • 



ds 

dK 

K 



- K 



IK 



(2-15) 
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and since H ■ q + Kp, then 



J - - K 



iH 

n. 

UL 



*p 

dq + K dp 
ds ds 



- - _Ke_ 
is 



S«-8, 



- - 0i 



s—s, 



1 \ 

The relation S R J « -Kp/(q* + Kpv is quite convenient since it only 

requires knowledge of the "open loop" transfer function and need not 
be in factored form. The evaluation of the pole sensitivity S^ J by 
the residue of the closed loop transfer function will call for the 
complete determination of the closed loop transfer function. 

II-4 - Slope of Root Locus 

Slope of the root locus is the tangent of the argument of s 



and since, 

s 8j - 

S K 



Then As^ • 



As 



1 



AK 

K 



8 



j 



K 



AK 

K 



j 

(2-16) 

(2-17) 



Since AK is a positive number, then 
K 

'As. 



± 



A' 



where A 6 * represents arg As^ 



*2 m 
dx 



So 

ta n ten ^ 



w 

ReS®j 



(2-18) 



(2-19) 



But from Equation (2 



w 






i (2- 


>8) we have, 


i .. . . 1 


1 




2 1 


- 2 1 




I 

m 




i?i * + ‘i J 


1 S»-8 






1 n 


1 




2 1 


- y l 




Re 


j-i 4 + p j 


i?i 8 + z i J 


8—8 



(2-20) 
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Input 



-e- 




Output 



Fig. 2-1: Unity Feedback System. 
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If we wfrite 



s + Pj " 

8 + z i « i ± e^ 0 i 



and since 



-z- 



sin 9 



1 



then 



m sin 9. 

t L 



n sin 9. 



dx 



>1 tl 



i-1 



m cos 9 



2 -j- 

j-1 7 j 



1 



n cos 9. 

2 ~r 

i-1 



Equation (2*21) could be written in simple form, 

. 2 A i 



dx 

dx 



sin 9. 

TT 1 



Z a C08 0, 

A i“ 77 

where - 1 if J ^ and 9^ are due to poles 

A^ - -1 if and 9^ ore due to zeros. 



and 



Then 



- magnitude of -8^ + 
9^ - argument of -s ^ + s^ 




( 2 - 21 ) 



( 2 - 22 ) 



(2-23) 



In some cases it is easier to compute dy/dx from the slope of ds/dK 
obtained by differentiation of - K - p/q as shown in example discussed 
at the end of this chapter. 
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II-5 



8 

Graphical Evaluation of 



J 



The pole sensitivity at any point "A" corresponding to a particular 
value of s can be evaluated graphically by summation of angles and with* 
out evaluation of K. Choose a point C such that AC is perpendicular 
to the locus and evaluate the phasor AC ■ AS Evaluate the value of 
the phase of K at C by summing angles ©^ and let it be "Tl + Aol This 
procedure is based upon the Independence of the derivative of an 
analytic function on the direction of differentiation. Hence AK is 
chosen to be A'C' which is mapped into AC. 

Then: K - j(A*)K 



*1 

As - As « As ■ S J (2-24) 

i** i 4K 45 

JK K 



Therefore knowing the pahsor As and the angle A* one can evaluate 
II-6 - Curvature, Pole Sensitvity and Their Relation 

The curvature at any point on the root locus is proved to be 
(Appendix II), 



1 



Im<^ K >ss 

Imj(l*K). 



Im<^« 
Re(jUK). 



(2-25) 



Z 2 

At first glance it seems that the curvature and sensitivity are somehow 
related. An example, discussed at the end of this chapter, will reveal, 
however that the curvature in that example is constant while the pole 
sensitivity is equal to cot <|> and consequently is a continuous variable. 
The independent nature of the curvature and pole sensitivity can be 
understood by noting Equation (2-25) and the definition of sensitivity. 
Pole sensitivity depends on first derivatives, while the curvature 
depends on the second derivatives, and in general the two derivatives 
are independent. As suggested by Ur we can think of curvature as 
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"kinematic" property and the pole sensitivity as a "dynamic 
property." 



II-7 - Open Loop Pole and Zero Sensitivity (Figure 2-5) 

Let G(s) be the forward transfer function and L be the loop 
transfer function. Then L » GF and the system characteristic equa- 
tion is 



1 + L » 0 



and if -Sj is the system root, then 



1 + L(s) 



8*®8 . 



where s » -s^ represents the point on the root 
gain constant of L(s) and Z^, P^. its open loop 
we can write - 



locus. If K is the 
seros and poles, then 



L - L(s,K,Z i ,P i ) 



and taking the total differential of L; 

n 

2 bL_ 



dL - bL Att + bL Av + 

u * Sit dK 



i-1 iZ, 



dZ, 



m 

+ 2 U 

i«l bP, 



dPj 



(2-26) 



On the root locus 1 + L » 0 or L(s) » -1 and therefore L(s) ■ constant, 
and the total differential dL is zero for s » -s^. Then letting 
dL ■ 0 and s » -s^. Equation (2-26) gives. 



0 



Jk 

bs 



8*-S 



dSj + 



bL 

iK 



dK + 



jj® 



j 




8“ -8 





Rearranging we get 



ds, 




dK 



s»-s . 






(2-27) 



But itself is a function of K, Z^ and P^, so 
8^ » 8j(K, Z v P t ) 



( 2 - 28 ) 
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Taking the total differential of we get 



ds . 

J 



K i!i ^ + 2 — i dz i + 2 S dP 

Uk K + i \7,. 1 + i SpT 1 



(2-29) 



Equation (2-29) suggests that ds^ be written as 



ds 



S S 8 

fi + Z h[ " Z 1 + 2 S pJ dP l 



(2-30) 



Equations (2-29) and (2-30) thus define the sensitivity of root s 



j 



to gain K, as 



s! J 2 > 8 jl 



K 



(2-31) 



which is the same as Equation (2-13) and the sensitivity of root e 



to open loop zero Z^ as 



5 j D bs. 



bZ, 



j 

(2-32) 



and the sensitivity of root s^ to open loop pole as 

< J 2 ^ 

S . s , s . 

II-8 - Relationship between S„ , and S J 

K u F 

Comparing Equations (2-27) and (2-30) yields: 



(2-33) 



K 



k f hAk \ 

\ ^LAS / 8 ■ -s 



j 



^ 8 • -8^ 

( ^ 

\~Tl7Ts / 8 ■ -s. 



Tl7Ss 

^L/fcP 



(2-34) 

(2-35) 

(2-36) 



But 



TT (s + Z.) 
L « K 1 1 



m 



TT (s + p ) 

l 1 



(2-37) 
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Then 



TT (s + z.) 



(ik) . _i 1 

l fcK ) s m ~B , m 

v 1 7T<. + p t ) 



(k )' - ->j 



(2-38) 



S « -8 , 



But since 1 + L(s) * 0, 

1 8 ■ -8 



then L(s) 



I 8 ■ -8, 



Then (2-34) become® 



Similarly 



J 

1, and Equation (2-38) yields to 



(i). 



1 

K 



• • V(s) . - 






n-1 

kIT (s + z ± ) 

IT (s + p ) 
l 1 



8 ** »S. 



(2-39) 



(2-40) 



8 + Z 



i I 8 " ”8, 



* 8 j + Z i 



(2-41) 



Then, 



OR 



UAs 



s rJ 

- 8j + z t 



S • °8 . 



•l/(- »j + z t ) 



(bL/bS) 



S ■ -8 



j 



(2-42) 



Similarly 



Sj - P A 



(2-43) 



Equations (2-42) and (2-43) dictate that the sensitivities to all 

8 i 

singularities are directly proportional to S^ and inversely pro- 
portional to the distance between -s^ and the singularity involved 
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(in this case s » - Z ^ and S - -P are the signularities) Thus what- 

8 1 8 1 8 1 
ever properties are found for may be extended to or Sp 

A particular case of Equation (2-43) is for P q * 0, i.e. the 

root sensitivity to the pole at origin of the S-plan. 

8 . 



S 



j 



j 



(2-44) 



‘•j’ 



S 1 8 1 

Then S J is proportional to S„ , the constant of proportionality being 

o 

1 . a complex quantity. 

II-9 - Relationship between Root Sensitivity of Sj and Residue 

In Equation (2-14), we proved that in case of unity feedback, 

s, 



j 



K 



-Q 



j 



(2-45) 



where Qj is the residue of T(s) 



KP(S) 



q(s) + Kp(s) 



at point -s 



j' 



In this section it will be shown that if -8j is a single root, then 



F< -*j) 



(2-46) 



where F(-s^) is the feedback transfer function F evaluated at s * -s j 
and Qj is the residue of the system transfer function at (-Sj). For 
unity feedback and simple root F(-Sj) » -1 and Equation (2-46) yields 
to Equation (2-45). The overall transfer function T(s) of the system 
shown in Figure (2-5) is 



T(s) 






1 + F(s)G(s) 
But L(s) *» F(s)G(s), then 






F(s)[l 4- F(s)G(s) ] 



(2-47) 



T(s) 



L(s) 



F(s)fl + L(s)) 



The residue of T(s) at -s ^ is 



Qj * (s 4- Sj)T(s) 



8 ■ 



(s + s j)L(s) 
F(s)f l 4- L(s)l 



(2-48) 
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Then, let the rightmost expression be denoted as R^, i.e. by 
definition. 



Then 



V'V - Qj 

RjF(l + L) ■ (s + s j)L 



(2-49) 



Taking the derivative with respect to s of both sides, 

[1 + L) Rj + F Rj + F t l + L) - L * (s + Sj ) ^ 

58 

At s ■ “8. i.e. at a point on the root locus, L * -1 the above 

equation reduces to 
IL 



FR 



OR 



j Is 



y-y 



8 « -S 



j 

1/F(=a 



> r*i 

L * 8 J 8 " -Sj 



Comparing Equations (2-50) and (2=40) we get 

4 J - F <*y«j 

Now turning to the case of N order roof at -s^ we can expand 



(2-50) 



(2-51) 



T(s) 



F(1 + L) 
into partial fraction as 



T(s) 



fl3L 



S + 8 



j 



where 






(N - K)! 



+ 8,)^ + 






y 



+ terms from other roots 



(s + Sj) 



N 



vN-K 



N-K 



(s + Sj)^ 
F(1 + L) 



(2-52) 



(2-53) 



Again defining the quantity inside the small bracket as R , then 
repeating the operations as before, one obtains after repeated 
differentiation 
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(2-54) 



jN 



- N! 

F(-S ) 

J N 



8 ■ -8 



and from Appendix III 

s *i - 

S K 

Theorem 



(-1) N N.' 



(2-55) 



[ 5 ] 



8 « -S 



When -8j is a single order system root, the sum of the 
sensitivities of s^ to open loop poles and zeros is equal to unity, 



i.e. 



Z £ + 2 s J 



i z i 



i p i 



(2-56) 



This is easily seen by referring to the construction of root loci. 

If all open loop zeros and poles are displaced by the same amount 
T, then all closed loop roots are displaced by the same amount ^ , 
that is if dZ^ • dP^ for all i, then dSj for all i. There- 
fore if in Equation (2-29) we put ds^ » dZ^ ■ dP^ ■ 5 and dK « 0 we get 




- 1 



A more rigorous proof is as follows. Rewriting equations (2-13) 
and (2-42) and (2-43) we get 
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Then n a 

I s J 



Ff 


1 , 


m i 

1 




[i-i z i 




°ri 

0 

00 

•j H 




m s. 

+ 2 s J 

i»l i 


n 

- 2 
i«i 


£ 

s J 

- i V 


s *i 

S K 


Z «= s ^ ^ 

i j t-i 


8 j * P i 




8 j 

* V 


n . 

1 , 


m 




^ Z ® s 
i»l Z i j 


^ s 
t-i 8 j 




- i 







The above relation Is only valid for simple poles. 

11-10 - Graphical Method for Determining Open Loop Pole and 

Zero Sensitivities 

Equations (2-42) and (2-43) suggest that to singularities which 
are close to system root (-s^), ®j more sensitive and for singular- 
ities which are far away, s ( is less sensitive until it becomes 
insensitive to singularities at infinity. This concerns the magni- 
tude of sensitivities. But sensitivities are vector quantities, 
since the change of a parameter may move the roots in different 
directions. It is then helpful to make use of Equation (2-56) toget- 
her with Equations (2-42) and (2-43). 

Referring to Equations (2-42) and (2-43) where Z and P indicate 
open loop singularities and -s , is the system root in question, we 
can draw a vector from -s^ toward each and every pole, and away from 
each and every zero. The length of each vector is inversely propor- 
tional to the distance from -s^ to the singularity concerned. Then 
construct the sum U of all these vectors, which is a vector itself, 
(Figure 2-5). If U is taken as a unity vector in magnitude and phase, 



then the other vectors measure the sensitivity of to each 
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S-Plane 





Fig. 2-5: (a) Construction of vector diagram 



(b) Vector addition to obtain U. 



30 




Fig. 2-6: A Basic Feedback System. 
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singularity respectively. So we regard it as 



U - l /o - 1 



J° 



(2-57) 



and is named by Thaler and Rung as "unit-sensitivity" vector. This 

8 i 

unity scale does not apply to sensitivity to gain, S^ J . Actually 
this technique is quite similar to the graphical method of evaluating 
the residues for simple singularities, and by putting U « 1 /o_ we 

establish the necessary reference scale. Phase of sensitivity vector 
is measured as positive in the clockwise sense starting from the U 
vector. This seemingly arbitrary sign convention in fact comes 
from Equations (2-42) and (2-43) which are the basis of this vector 
diagram. 



s s. s. 

S*J - V - V ■ V (2-58) 

1 v P t <-v • '-)> x 



The denominator is the vector from root (-8^) toward pole (-P^) 88 
shown in Figure 2-7. The phase^ relationship of the above equation is 

'j 



/**! ■ / S K J ■ Za 



(2-59) 



Relation (2-59) is true no matter what conventions are applied to the 

measurement of the angles. Since the U vector has been found to be 

8 8 8 

equal to 1 /o as far as S_ J and S, are concerned, phase of S_^ 

l Z 1 P i 

will be measured starting from U as zero phase. On the other hand 
s . 



/ s 8j / ~V 

/ K and j V 1 



are measured in the conventional way, i.e. 



starting from the positive real axis and counting positively 
counterclockwise. 

In Equation (2-59) the quantity / "'K does not depend 






on i, i.e. it is the same for all i's. Thus for each i, 

» constant 



xo biic c 

Zi 



A 
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meaning that the larger the value 
/ Sp l . Since / V 1 

A 



of / \ 



the smaller must be 



is measured conventionally (positive 



counterclockwise) then 
wise. 



Fj must be measured positively clock- 



s . 
J 



Next we prove that vector S J lies on U. From Equation (2-44) 

K. 



S , 

s J = 



K 



(s ,) 



and 



lA - ZlL + Zi 

Figure 2-7 shows the angles / ’j and J S ^ 

L — !° 



= ct , the latter 
s , 



being measured from U. Since as shown in Figure 2 



/ S j 

-7, / P o + 



Lh l 



gives the direction of U on the S-plane, it is thus establish- 



s s . 

ed that on the S-plane, S ^ always lies on U. Since S J indicates 

K. K. 

the direction in which the root moves when K varies, i.e. the direc- 
tion of the root locus, the above result can be stated as: "At any 

point on the root locus, the U vector is tangent to the root locus." 

Also Equation (2-59) reduces to 

s 



LL ■ 'Zi 



(2-60) 



11-11 - Locus of U on the S-Plane 



Consider the plant to be compensated. 



G = 



K(s + *1) 



s ( s + P 1 )(s + P 2 ) 



(2-61) 



and the desired dominant-root location -s^ as indicated in Figure 

2-8. A spirule measurement shows that an additional phase of + <{> 

is needed at location -s . . Thus a lead network with a zero at Z and 

J 

a pole at P is needed. The question is, how does the tip of the U 
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vector move on the S -plane, when Z and P take all possible values 
on the negative real axis? 

Now sensitivity vectors should be drawn for all the plants poles 

and zeros, ignoring for the moment, the need for compensation. All 

these vectors can then be added to obtain the vector sum QI, which 

may be called the "uncompensated" unity vector. Then the geometric 

locus of U is a circle with a center I and a radius of r = 1 . 

d sin (J) 

This circle is named the (U) circle. 

Proof : 

QI would be the unity vector without "P" and "Z" (Figure 2-8a). 

a , 

When P and Z are added so that PQZ = <j», a sensitivity vector must 

be assigned to each of these singularities, Qm = 1 is assigned to 

QP 

the pole and Qu = 1 is assigned to the zero. The vector QI then 

QZ 

increases by the vector quantities IM = Qm and MU = Qu. QU is the 
final U vector. Note that, in Figure 2-8a, IN and MU are two equal 
but opposite vectors. 

Figure 2-8b shows that when P moves along the real axis, since 

Qm = _1 , m moves on a radius of r = _1^, because of inversion of 

QP 2d 

real axis with Q as center and ratio of inversion one. Such a circle 

is referred to as the (m) circle. If we translate this (m) circle 

by QI, then we get (M) circle which is the locus of point M. So the 

locus of point M is the (M) circle, or radius R = _1 and with I as 

2d 

its upper most point. Figure 2-8c shows that: 

IU = MN = 2R sin MON 

2 



IU = 2R sin <J) = 1, sin <f) = constant (2-62) 

d 

Hence U moves on a circle centered at I and of radius r = 1^ sin (|. 

d 



So the following result can be stated: 
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u 




Point U. 



35 



"The geometric locus of the tip of the unity vector Is a circle, 

centered at 1 and of radius r ■ ^ 8in 4> where I is the tip of the 

d 

'uncompensated' unity vector, d the imaginary part of the dominant 
system root and (j) the phase shift to be Introduced by the compensation." 

Note that there are limits to the geometric locus of U on the 
(U) circle. One limit corresponds to the extreme case where P ■ <», 
the other is for Z * 0. Later it will be shown how to obtain these 
limit points on the (U) circle. 

11-12 - "Use of the U Locus" 

Starting with the uncompensated plant and a desired location for 
the dominant root, it is always possible to be able to draw the (U) 
circle and be able to obtain the two limit points on the U locus on 
the circle. 

Conversely after a point U has been chosen on the U locus, it 
is possible to perform the graphical construction in reverse order 
and thus drive the corresponding Z and P, i.e., obtain the compensation 
needed. The following will explain how to choose the U point on the 
U locus, which is the seence of this construction technique. 

In Figure 28a, note that OJ is perpendicular to IU where J is 
the midpoint of arc MN. P and Z can be obtained by drawing OJ per- 
pendicular to IU, cutting the (M) circle at J, drawing arcs JM ■ JN ■ 

<f>. P is then determined by drawing OP parallel to IM, and Z is det- 
ermined by drawing QZ parallel to IN. 

An example is given in Figure 2-9a, where a desired U point is 
shown (an arbitrary choice for illustrative purposes). OJ is drawn 
perpendicular to IU. Next, arc JM » arc JN - (j) is drawn. In Figure 
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t® 




u circle 
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U circle 



1 




(b) 




\ 



/ 



/ 



Fig. 2-10: Determination of the limit points of geometric 

locus of U. 
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2-9b P is obtained by drawing QP parallel to IM and Z is obtained by 
drawing QZ parallel to IN. 

Now the problem recurs of determining the limit points of the 
U locus, using the same construction procedure. The extreme right 
limit point corresponds to the case where one compensator's singu- 
larity (here, the zero) is at the origin. The extreme left limit 
point corresponds to the case where the other singularity is at - oo. 
Figure 2-10a shows how to obtain the right limit point, Draw- 

ing IN r parallel to QO, because Z is now at 0 and take, arc N r J r = 
Drawing IU^ perpendicular to OJ^, then results in U r 

Figure 2-10b shows how to obtain the left limit point, U^. 

Draw IM^ parallel to QP, which is horizontal since P is at - oo, and 

take arc = <j* - Drawing IU^ perpendicular to OJ^, then results 

in U^. The geometric locus of U is the part of (U) circle between 

U. and U . 

1 r 

11-13 - Design Techniques 

In the preceeding section it has been derived a method for 

finding P and Z, given the desired location of point U on its locus, 

i.e. the magnitude and phase of QU, the unity vector. This method 

is restated below in a step by step form. Figure 2-8a. 

Step 1 - Considering Q as if it already were a point on the root 

locus, draw the vector diagram at Q for the uncompensated 

system and obtain QI, the "uncompensated" unity vector. 

Step 2 - Draw the circle diagram, composed of the (M) circle, of 

radius R = _1 and whose uppermost point is I, and the (U) 
2d 

circle centered at I and of radius r = 1^ sin (j>. Fix the 

d 
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limits of the locus of U on the (U) circle. 

Step 3 - Given QU as the desired unity vector for the compensated 
system, draw OJ penperpendicular to IU which cuts (M) 
circle at J. Draw angles JON = <{). Then Z and P are deter- 
mined by drawing QZ parallel to IN and QP parallel to IM. 

In this section, it will be seen how the "U n vector (i.e. the 
location of U on its locus) is selected to satisfy a particular 
condition. 

(a) Design for Minimum Root-Sensitivity. 

Since the U-vector is the scale used to measure the individual 
sensitivity vectors, the larger the scale, the smaller the magnitude 
of the sensitivity measures. Thus one possibility is to design the 
compensation for maximum magnitude of the U vector, that is minimum 
root sensitivity to open loop singularities. Figure 2-8a shows 
that maximum magnitude of QU is obtained by placing point U near the 
lower most part of the (U) circle or more exactly on the extension 
of QI. If such point is not within the locus of U, then it can not 
be a location for U and one must select the lowest point which is on 
the locus. In Figure 2-9a, this is point U^. Thus QU^ is the selected 
unity vector, and with this given, one can proceed to the 3 step pro- 
cedure outlined in the beginning of this section. With such a design 
s . s . 

S ^ and S ^ are all minimum, for all i’s. 
i L i 

(b) Design for Constant Damping When K Varies. 

Another practical problem is to compensate a system in such a 
way that when gain K varies about its nominal value, the dynamic 
response of the system does not change. This calls for a constant 
i.e., a root locus that remains tangent to the radial line OQ at the 
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neighborhood of Q, Figure 2-8a. How can Z and P be found to obtain 
such a root locus? It is now shown that this can be done by merely 
selecting point U so that the unity vector QU goes through the origin, 
0, of the S-plane. In other words choose U so that Q, 0 and U are 
in line. If the locus of U does not permit such a choice, this means 
it is not possible to obtain a constant ^ about Q for the given system. 
One can then choose the best solution available, by taking the "U" 
location that is closest to a straight line with Q0. 

The above statement can be proven very simply if one recalls 
Equation (2-44) that s. 



s . 
! J 



K 



s . 
J 



(2-63) 



s . 



The specification here is to force S ^ to have the same direction as 

K 

0Q. Thus making s. move on a radial line when K varies. This means 

Z s . J 

S K J must be equal to the phase of Q0 which is also phase of s^. 
But from Equation (2-63) 



LA L ■ L 



s . 

S J 
S K 



b x 



= 0 



(2-64) 



g 

The phase of j must be zero, this means that the sensitivity vector 
s . Po 

S p J must be on the U-vector or conversely, the U vector must pass 



through P at the origin of the S-plane, 



s . 



A faster way to prove this is that we show that S ^ always lies 

^ K. 

on the U vector. In order to keep ^ constant, must be radial. 



thus U vector must be radial, 



11-14 - Limiting Behavior and Special Cases. 

The magnitudes of gain sensitivities can cover the entire range 
of values from minus to plus infinity. Yet, intuitive notion of 
"sensitivity" as a general concept in closed loop systems makes part 
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of this range unreasonable. One part of the problem is a direct 

consequence of the sensitivity definition, while another is associated 

with its first order approximation nature. A better understanding of 

both facets can be gained by an examination of limiting cases. 

In general, closed-loop poles depart from open-loop poles for 

low values of gain, and proceed to either open loop zeros or unbounded 

values as the open loop gain becomes very large. The gain sensitivity 

as given by Equation (2-13) could be written in the form 
s . 



K 



(2-65) 



2 — - — 
i=i z i - s j 



m+n 



+ X 



i=i s j 



As K approaches zero, the closed loop root ( -s j) approaches the open 

loop pole from which it derives, i.e. Sj-*P^. Then the term, the 

above relation reduces to 
s . 



K 



K-r 0 



( S j - P i ) 



+ 0 



( 2 - 66 ) 



Similarly as K becomes very large, n of the closed-loop poles approach 

open loop zeros. If the jth closed-loop pole is one of these, and 

it approaches the ith open loop zero so that 

s . > Z. 

J i 



Then, 




(2-67) 



Finally, m of the closed-loop poles have no zeros to go to, and hence 
become very large relative to P^ and Z^. The sensitivity for these 
poles is 



<5 




1 






i=l 



1 



s 



J 



S 1 ( 2 - 68 . 

m 
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Here m = (order of denominator) - (order of numerator). When the gain 
is sufficiently large for the open-loop zero db line to intersect 
the high frequency asymptote on Bode plot, the open loop transfer 
function is approximately 



L(s) 




where m = (order of denominator) - (order of numerator). Then for 
the points on root locus 



1 + L(s) 



0 



Is = -s . 

J 



or 



1 + 



K 



( - S j ) 



m 



or 8 = - y~- (2-69) 

Thus, the sensitivity of the unbounded pole from (2-68) will be 
s 



K J s j >> z i p i 



m 



¥2L 



(2-70) 



m 



Equation (2-70) indicates that sensitivity increases as the m root 
of K, as gain is increased, although for the finite gains the sens- 
itivity is always finite. Referring to Figure 2-11, the above special 
cases can be simplfied as: s 



For poles going towards "O.L. zeros" 



For poles going towards infinity 



when K -*■ oo 
when K 0 



s . 

S J = S 
S K 



K 



0 



4 j ^° 



in general 



m 



S KJ = 



m 



y~~K when L ® K 



m 



m 



Another circumstance in which the sensitivity can become very 



large is revealed by Equation (2-45), and that happens when (the 
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Fig. 2-11: 



Poles going Cowards open 
loop zeros. 




Fig. 2-12: Poles going Cowards infiniCy. 
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residue evaluated at pole -s^) is not finite. This is to be expected 
since the sensitivity factors defined thus far have not considered 
multiple-order, closed-loop roofs. As far as the gain is finite, 
an infinite gain sensitivity always indicates multiple-order closed- 
loop poles. 

A special situation of considerable interest can occur when a 
closed loop root lies between an open-loop pole and zero which are 
much closer to each other than to all other open-loop poles and zeros. 
This is the so called dipole case as in Figure 2-12. The sensiti- 
vity for the bounded closed-loop pole will be, approximately 






<W<W 

Z 1 - p l 

(2-71) 



Differentiating equation (2-71) to find values of s. which make, 

J o 



sensitivity maximum, reveals that maximum value of will occur 

when s , _____ 

J 2 

8 , 



Z + P 4 

i i for which S r J becomes 

2 



j 



max 



1 (Z. - P*) 
4 1 1 



(2-72) 



11-15 - Sensitivity at Irregular Points 

Expanding the gain K * K(s) in series to include the higher 
order terms we have 

K - K + dK As + 1 d\ As 2 + — — (2-73) 

° d » 1 ds 2 

and then assume the first n°l derivatives vanish* This means then, 

AK-K-K » 1 d°K As n 

^ fl* «l 
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or 




A 8 ■ n Ini Ak 

/ d°K 



d°iC / d°K V^K (2-74) 



Since As 



V2K 



-4 

is finite it is suggested to be substituted for S^ J 



at irregular points. 

For anyAK> 0, ds has n values separated by 2TT which represents 

n 

the outgoing branches in Figure 2-13. If AK < 0, n branches with the 
same separation as before are again obtained, but they are displaced 
from the previous set by Tf/n, that is the two set Intersect. This 
could be stated as: The incoming branches, which at their junction 

represent an nth order closed-loop pole, are evenly spaced and sep- 
arated from each other by 277/n. The outgoing branches are also 
separated from each other by 277/n and are midway between the incoming 
branches . 

11-16 - Sensitivity Functions for Alternate Transfer Function Forms 

The transfer function could be written either in root locus or 
Bode form. For Figure 2-6, open-loop transfer function can be written 



as 





+ + a n 




n 




(2-75a) 



It <® + V 



i-1 




(2-75b) 
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The Equations (2-75a) and (2=75b) represent the root locus and Bode 
fora respectively. 

The equations for the gain sensitivity are the same whether open- 
loop transfer function is in root locus or Bode fora, i.e. 

s , 



j 



K 



S j 
9 



(2-76) 



This is because in comparing (2-75a) and (2«75b) it is obvious that 

^ ® K x constant 



and then. 



K 




It is however necessary to modify the open loop pole and zero sensi- 
tivities for terms which are written in Bode fora. For poles and 
zeros which appear in open loop transfer function in Bode fora, i.e. 

, the sensitivities are 






or 



+ 1 



’j 



3 1 S K 



j 



V Z i “ S j> 



and 



8 1 V 



V 8 j - V 



(2-77) 



(2-78) 



Frequently open-loop zeros and poles will occur as complex 
conjugate pairs, and variations in the system will change both zeros 
or poles. For example, consider a complex pair of zeros, Z^ and Zj> 
which are defined by their frequency, is), and damping ratio ^ , i.e. 

Z, 



ScO + jW\/l - 

- jtt/l - ^ 2 



(2-79) 
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For this situation it is useful to define frequency and damping 



ratio sensitivities as 



& - 


& 


- !!i 

iZ. 


dZi 

dco 


♦ t!i 

bz. 


dZ 2 

d 


S 8 j - 




b z i 


dz i 


*** 


dZ 2 



(2-80) 



Where ^ 8 j and ^ 8 j are given according to (2-42) or (2-43) (In case of 



iZl iZ2 
poles). Components as 



dZ„ 



“i «i «2 -2 

TET. <Tf- 75" ■ * n<l rf can ' a,lly b ' * v * 1 - 



uated from (2-79). Then when the complex pair of zeros appear in open- 

2 2 

loop transfer functions in root locus form, i.e. (s +2^4fe+»), the fre- 
quency and damping xatio sensitivities for a complex pair of poles or 

,s, 



zeros are 8 j _ -2,(6) 7 s* 

s] - 2^ i +4? 



j 






2 - 2$6>s, + 



(2-81) 



J g j 

Where the upper sign is to be used for zeros and the lower one for poles. 
If the term is written in Bode form, 



sitivltles are 



6 i + 2s j 

v - V 1 



lur zeri ana tne lower 

», i.e. I* + 2$s j| 

a - W J 



the sen* 



’j 



4* 2s i“ s k 



j 



s 2 - 25^8 J + W 2 



(2-82) 



Sj - 2^8 j + 

For some cases it may be more convenient to define a complex pair in 
terms of their real and imaginary parts, i.e. 



Z^ ■ a + jb 



a - jb 



2 2 2 

In this case with the terms in root locus form [s + 2as + (a + b )), 



the sensitivities are 



8 



j 



a 



da 



U 



\s 






dZ, 



da + ^Z 



dZ 



_2 

da 



Z 1 ’ 8 j 



( 1 ) + 



K 



(1) 



r * + i i 

(a + jb - 8^) (a - jb - s^) 
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a 



2(a - s ^ ) S, 



j 



2 2 
2aSj + a + b 



’j 



and in a more general form as in case of Equation (2-81) it could be 
written as s , 



’j 



a 

8 , 



"b a 

it 



+2(a - 8^) 



j 



2 2 
- 2as . + a + b* 

+ 2b S® j 



’j 



(2-83) 



2 2 

2aSj + a + b 



where the upper sign is to be used for seros and the lower one 

/ 2 2 2 

for poles. If the term is in Bode form l(s/(a + b )] + 
(2as/(a 2 + b 2 )] + 1^ , the sensitivities 



are 




(2-84) 



11-17 - Example 

In Figure 2-14, suppose the gain K, is given as 

1 



K 



2 

8+1 



Find radius of curvature, and pole sensitivity. 



He can write 



f. - + 7 



d 2 K 

ds 2 

d 3 K 

ds 3 



2 

3 



_6 

4 



Double points occur atdK^ 0 or S ■ + 1. Substituting S 



i 



ds 
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Fig. 2-13: Outgoing Branches (n ■ 2) 



y 




we get 



dK 

ds 



— a -1 + 



= 1 + 



Jff 

°2J0 



2je°”^ e^ ° e‘ 



2j 

-2je ^ sin 0 
2 sin 0 e“ j( ^ + 2.) 



On circular portion of the root locus. 



T dK 
. In t- 

dy ds 

dx _ dK 

Re 3— 
ds 



The pole sensitivity is 
£ ■ *£ ■ 



cot 0 ® tan (0 - T ) 

2 



+ 1 




1 

8 



On the circular portion 

8 * J0 + e° J ^ 



j 



$4 



+ e r , JW ^ /i 

J 0 " -j jr “ ” j ® cot r 



e ' = e 

The pole sensitivity at s ■ -1 (a double point) is infinite and hence 
the suggested substitute, i.e. 
then 



A s is evaluated. Here n ® 2 and 

V*K 



A s 

Vak 




e + 1 evaluated at s « -1 
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Ill 



SENSITIVITY AND LARGE PARAMETER VARIATIONS 
III-l - Discussion 

This section deals with the most practical aspect o£ the sen- 
sitivity and is quite utilized in design consideration of the feed- 
back control systems. It discusses sensitity in two particular cases 
one with leakage transmission, and the other case when leakage trans- 
mission is zero. It relates the sensitivity of the system to return 
difference and null return difference. 

In previous sections it has been discussed that in most cases 
we are facing the variation in the values of an element in an 
engineering system. This element could be a vacuum tube, transistor 
or sometimes could just be a passive element. It is possible by 
means of feedback around the troublesome element, to achieve the 
desired reduction in sensitivity. Let the troublesome parameter be 
k, as in Appendix 1. We relate k to the controlled source S and 



control variable C, l.e. S • kC and, 
T 



kt .t 

t + c .i-° 8 _ 

oi 1 - kt 



(3-1) 



C8 



with the fundamental signal flow graph as shown in Appendix 1. 
We defined the classical or system sensitivity as. 



AM_ 



>T/T 



(3-2) 



'k i£«.k dk/k 

T 

so is the relative change in T divided by the relative change in 
k for infinitesimally small changes In T and k only. 

Case 1; 

Consider the case v;here t Q j« T, i.e. leakage transmission is 

negligible. T! is is a situation in many feedback systems. By straight 
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forward differentiation of Equation (3~1) we get 



„T \SU,T i T/T k 

s k s TOT * TkTiT t TiT 



kt ,t 

ci 03 

1 - kt 



cs 



t ,t (1 ~ kt ) + t (kt ,t ) 
ci os cs cs ci os 

( 1 - kt ) 2 

CS 



kt cl t o e 



t .t 
ci 08 



(1 = kt )' 

' C8 



and from that 

s k - T^Tt — * -r- for *oi ’ 0 < 3 ‘ 3 > 

cs k 

where F. is the return difference' for reference k. It is clear now 
k 

why It is insisted that all t^’s i* 1 the fundamental feedback equation 
should be independent of k. To have a small sensitivity then | 
should be large, i,e. (-kt ) which is usually called the loop gain 

CS 

should be a large number. 

T 

Since t and possibly k is a function of frequency, so is S. , 

CS K 

T 

and one must consider the frequency range over which a small is 
desired. It is interesting to note that as the magnitude of the loop 
transmission determines the sensitivity, its phase also determines 
the stability of the system. 

Case 2: 



When the leakage transmission t^ 4 0, then by straight forward 
differentiation as in Case 1 we get 



f 




Fig. 3-1: Signal Flow Graph for Example (1). 





Fig. 3-3: Signal Flow Graph for Example (3). 
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This could be written in the form 



S. 



(1 - kt ) 
cs 



kt .t 
ci os 



kt ) 
cs 



and since from Equation (3=1) 
T - t 



kt .t 
ci os 



oi 



kt 



cs 



then 



kt 



cs 



[■-*] 

] 



t 

— 2i_ 
T 



(3-4) 



S, can also be written in terms of return difference and null return 
k 

difference. From Appendix 1, Equation (A3-3), we have 




(3-5) 



where F^ is the return difference with reference to k and F£ is the 
null return difference also with reference to k. Equation (3-4) 
suggests a practical mean to measure the sensitivity. It is noted 
that in case t^j » 0, Equation (3-4) reduces to Equation (3-2). The 
applications of Equation (3-4) are illustrated by a number of simple 
examples . 

Example - 1 

In the single-loop system Figure 3-1, if the sensitivity desired 
is with respect to the forward transmission, the leakage transmission 
is zero, because with K « 0 there is no transmission from input to 
output. The sensitivity in this case is the reciprocal of the return 
difference, i.e. 



(3-6) 




1 



1 



F, 



K 



’ T'-'/Tk 



Example - 2 

In the same single-loop system, if the sensitivity of interest 
is with respect to the feedback transmission , ft , the leakage trans- 
mission will then be K, and Equation (3-5) yields 



In this example it is noted that when ^ ■ 0, i.e., there is no feed- 
back around k, then the leakage transmission t Qi -lxKxl-K. 
Comparison of Equations (3-6) and (3-7) demonstrates the well known 
fact that a value of K large with respect to unity results insensitive 
to changes in K, but with the closed-loop gain essentially equal to 

1 (when K is very large) the sensitivity of T with respect to P 

f , o 

approximates the value of -1 » 1 /18 . 

Example - 3 

Figures 3-2 and 3-3 represent alternative designs for a feed- 
back amplifier. Three amplifier stages are available, each with a 
gain a, and the overall gain to be realized is specified (less than 

3 

a ). Negative feedback is to be introduced to reduce the sensitivity 
of the overall gain to variations in supply voltage. To a first 
approximation this voltage variation can be considered equivalent to 
a variation in each value of a. The question to be answered is which 
configuration gives a lower value of sensitivity of the overall gain 
with respect to a? For Figure 3-2, the overall gain T^ is 



The sensitivity of T^ with respect to a is three times the sensitivity 




(3-7) 



T 



1 
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of an Individual stage with feedback, or 



T 

S 

a 



1 



3 



£.a 



(3-8) 



For the second configuration, with a single overall feedback path, 
the overall gain T. is 3 

4 S 



and 




(3-9) 



The two sensitivites (3-8) and (3-9) are to be compared on the basis 

2 

of the equality of and T , This equality establishes the relation 

(1 - /^a) 3 si- ^a 3 (3-10) 

substituting (3-10) into (3-9) we get 

T 

s 2 - 3 —^ 

® (1 - ^ a) 3 

Thus, the single overall feedback path results in a system with less 
sensitivity, to changes in a. 



Ill -2 = Large Parameter Changes. 

A serious shortcoming of Equations (3-1), (3-2), and (3-3) is 
that they apply only for infinitisimally small changes of k. There 
is therefore uncertainty as to their applicability for moderate or 
large changes in k. It is therefore found useful to use a new 
definition, defined as follows; 

Let T , k represent the nominal or original design value of 
o o 

the system transfer function and of the element under consideration 
respectively, and let T^, k^ be the corresponding (final) values 
at the new value of k. Thus, 
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T + &T 

O 



k, - k + Ak 
£ o 



(3-11) 



The new sensitivity function is defined as 
AT/T„ 



T D 
S k " 



Ak/k. 



(3-12) 



Case 1 : 

When it is assumed that t 



f 

,T 



k t .t 
o ci os 

1 - k t 

O C8 

k.t . t 
f ci 08 

1 - k t 
f cs 



AT/T 1 

Ak/k, 



oi 



T. - T 
f o 



k. - k 
f o 



0, then from Equation (3-1) , 

(3-13) 

(3-14) 

(3-15) 



Substituting (3-13) and (3-14) into (3-15) we get 

, T 



k 



V 1 ~ k f t C8 ) 

V 1 * Vcs> 



k. - k 
f o 



1 - k t 
o cs 



1 + L 



(3-16) 



where F, is the initial return difference and L is also the initial 
k o ° 

loop transmission. It is thus shown that when the leakage transmis- 
sion t is zero, the classical and new definition have the same 
oi 

value, although they are defined differently. With the new definition 

there is no uncertainty as to the effect of large parameter variation. 

Equation (3-16) is exact only if t is zero, but it may be used 

oi 

whenever T^. Fortunately t is zero for a very large class 

of feedback amplifiers. From Equations (3-11) and (3-16) we have 
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T 



Therefore 



and 



T 

_c 

T. 



£ 

JL 

k £ 
o cs 



AT 

T , 



S 



T Ak 

k k. 



“ - 1 - (1 



k f / 

W i ) 

k. M 1 = k t ) 

f \ OC8 ' 

(*f ' k ®‘«) ( 1 - Vc. ) 



£ 
+ L 



1 + L 



(3-1?) 



where L 



»k t 
o cs 



is the origins! or the nominal loop transmission 
and Equation (3=17) relates the original and final values of the 
closed loop transfer function which is of very practical design 
application. 

Case 2 

In previous cases we assumed t . was negligibly small in com- 

01 

parison with T and T_ over the range of variation of k. This is not 

O I 

always true. Suppose in a system whose signal flow graph represent- 
ation is that shown in Figure 3-4 parameter B varies substantially, 
and feedback is to be used to reduce the system sensitivity to the 
variations in D. If we choose B ■ k of the fundamental feedback 

equation, then t . may not be small. However the problem can be 
oi 

handled by method in Case 1, if we simply let k represent the entire 

element, i.e, k . a(B + ED) 

1 - AID 




Fig. 3-4: A Signal Flow Graph with Leakage Transmission. 




Fig. 3-5: Basic Feedback System. 



R E 


G(s) 




> 


i 













C 



Fig. 3-6: Basic Unity Feedback System. 
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The element k is then incorporated into the overall system feedback 
configuration . The problem then reduces to Case 1. 

In the above manner the majority of feedback problems can be put 
into zero leakage category. For those situations where a parallel 
path from the overall system Input nod® to the output node (from I 
to 0) exists, the following results are obtained: 

Since T 



t + kt .t /(I » kt ) 
oi ci osf cs 



using Equation (3-11) and notation (3-12) we get 



A T/T, 

j 

4k/k. 



(1 - k t ) 
v o cs 7 



( 



oi 



or 






oi 



Equation (3-18) can also be written as 

1 



S< T * t ol ) 



1 + L 



(3-18) 



(3-19) 



and from (3-19) by using Equation (3-17) w® get 



oi 



oi 



c /k. ) + L 
o f ' o 

1 + L 



or 



(1 + L ) 
o 



(VV" 



VV 



L + k /kf 

o o 



and 



At 

T 



T 



(1 ' ‘ol'V < l ' V k f> 

\ + <v k £ > 



0 - 20 ) 



III-3 - Further Investigation into the Meaning of Sensitivity Function, 
In prev ious discussions we defined the system sensitivity as 



S„ 



1 + L 



(3-21) 
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and T Is the 



where L is defined as the open loop transfer function. 



closed loop transfer function. Let us now investigate the value 

dT/T 

dL/L 

Referring to Figure 3-5 we have 



G 1 GH 

n — ■ i 





1 + 


GH 


H 1 + GH 


But 


L - GH 


and 


dL - HdG 




dT/T « L 


dT 


■ L 


dT 




dL/L T 


dL 


HT 


dH 



1 - 1 
1 + GH 1 + L 



Now by definition dT/T D 

dL/L " 

to open-loop perturbations, 
we get 

Also from relation 



_C i.e., sensitivity of the closed-loop 
s 0 

and comparing Equations (3-21) and (3-22) 




0-23) 



we can think of the sensitivity functions as the error transfer func- 
tion for a unity feedback system with forward elements GH, because 

R - C - E 



or 



RGH 

“ 1 + GH 



E 



then 



E 1 1 

R 1 + GH " 1 + L 



In control system, the error transfer function is always desired 
small. Ideally zero, and we may therefore observe that a small sensi- 
tivity function is desirable. 

Since the sensitivity function is dependent upon frequency and 
is a complex number for any arbitrary frequency we must further inquire 
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into its meaning. 



L b 



Let 

then 




o 




1 




or 



or 




(1 + Pcos 

1 + fees 0 
1 + 2 Pcos p + 



1 

0)+j( H sin fT 
, r sin P 

® J 1 1 r r ' r " " ^ 

i + 2Tco8 p + r 

(3-24) 



i 



)/i + 2Pcos 0 + P^ 




T sin 0 

i + r cos ~p 



(3-25) 



It is desirable to interpret both of the above equations from 
a physical view point. We then proceed in this direction. 

A - A Synthetic Interpretation of S 

— — ■ — I ' O 

We define the following four quantities: 



.M 

’M 

A 

A 

M 

’A 



s A 



D 


d | T 1 / 1 T | 


per unit increment 


in It| 


(3-26) 




d|L|/|L| 


per unit increment 


in | L I 


D 

IS 


d Ll 


» increment in - 




(3-27) 




d / L 


increment in / L 






D 


dtT|/!Tl 


per unit increment 


in I T ]_ 


(3-28) 


S 


d j _ L 


increment in ^ 


f L 


D 

S3 


d 


* increment in 


(3-29) 




diL|/lL| 


per unit increment 


in | L | 





Where | & I and L £ indicate the magnitude and phase angle of the 



complex number a, respectively. It is necessary to compute j T | 



and /T in terms of | Lj and ZL. in order to evaluate these functions 

i - r. J <* 

Therefore | L | « P and /L * 0 (3-30) 



and 



T 



(3-31) 



h i + rj0 



let H - |H|e J ^- 

therefore 

T - 



and G ■ | G le^^- 



(3-32) 



P 



^ y/l + 2 P coa 0 + r 2 



e j[ - tan" 1 P aln j - L m] 

i +rcos o —* 



(3-33) 

We are now In a position to compute the functions defined by Equations 
(3-26) to (3-29). Note that 

n - I Hi • \G| dP -|H|»d(G| (3-34) 

j> - /H + /G d0-d/G_ (3-35) 

From Equation (3-26) ve have 
M dlTl/lTl 



M 



d I L | / 1 L \ 

“ * Hl v/l + 2 Ceos 0 +P 2 



t Hl\/l + 2 P cos 0 +P 2 . 



dP 



(3-36) 



But | H | is not a function of the variable part of P then 

1 



s|J ■ ^1 + 2 T cos 0 +P 2 

M 



r ... . c — - i 

Ln/i + 2 P cos \ + P 2 J 
dP 



1 + P CQ8 



1 + 2 T cos 0 + F 2 
By similar technique we find 

„A 1 + P cos d 



(3-37) 



S A “ 



4 - 

e A 

5 • 

M 



i + 2 r cos ^ + p‘ 



p sin 0 



1 + 2Tco 8 0 + r 2 



- P sin ^ 



1 + 2 Pcos + P 2 



(3-38) 

(3-39) 

(3-40) 
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From above results it is shown that; 



-M 

S M 



M 

[•:] • < 



(3-42) 



(3-42) 



S M' and are all real numeric® which sr© desired small, 



ideally aero. To 



utility of th® sensitivity £unction s 



the mapping of th© loci of constant ^ and 3^ on th© polar plane 

(hence all four) will b@ 

M 

B - Mapping of 3^ Loci onto 



We require 



1 4» Pcos 4 
1 Hh 2 P CO® $ + P 



M 



(3-43) 

Where M is a constant. Changing to cartesian coordinates we get 



1 + x 



2 

1 + 2s + x + y 



M 



(3-44) 



or 



[* + » - 4 r>] 



2 2 
+ y 



(• 



■U* 



Thus the loci of constant 



*5 



are a 



m/ 

of circles with center 



(3-45) 



coordinates 

and 



x 



1 

2M 



- 1 and y s o 
radius * 



1 

m 



(3-46) 

(3-47) 



„A 



C - Mapping of S" Loci onto th® Polar Plane. 
M 



W® require 



P sin 6 



1 + 2 P cos 0 + P 

where Q is a constant. Changing t© cartesian coordinates we get 



(3-48) 



-y .. 

' 2 l 

+ 2x + x + y 



or 



+ iy 



+ i > 

2Q 



2 






(3-49) 

(3-50) 
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Thus the loci of constant S* are a family of circles vith center 

coordinates x ■ - 1, y ■ - 1 (3-51) 

2Q 



and radius ■ 1 (3-52) 

2Q 

Both families of circles derived for and are plotted in Figure 
3-7 and it is observed that the region of poorest sensitivity is 
near the (-1, 0) point, the sensitivity improves in every radial 
direction from that point. 

D - Happing of | j Loci onto the Polar Plane. 



From Equation (3-25) we require 
1 



\/l + 2 Pcos / + r 2 



- M 



(3-53) 



where M is a constant. In cartesian coordinate it is 

2 



2 2 
(x + 1 r + y 



'(-f) 



(3-54) 



The loci of constant j jare a family of circles with center 
coordinates, x ■ -1, y ■ 0 (3-55) 



and 



radius ■ 1 

M 



E - Mapping of 7 o Loci onto the Polar Plane. 

From Equation (3-25) we require that 

-1 P sin ft m c , 



tan 



1 + p COS J} 



(3-57) 



or in cartesian coordinate 



y ■ tan C' (1 + x) ■ C(1 + x) (3-58) 

The loci of constant phase angle of the sensitivity function are 
therefore straight lines through the point (-1, 0) with slope C. 

This family of lines comprises the orthogonal trajectories to the loci 
of constant j J . 
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Fig. 3-7: Loci of Constant 



F - Result of the Interpretation . 

For any point on the polar plane 

2 



(<y *(<)■• i < 



(3-59) 



where each term is evaluated at the point in question. Thus the 

M A 

magnitude sensitivity function is a vector measure of and at 

every point of the polar plane. 

The only extent to which the phase angle loci will be Of use 

comes from the consideration of the phase angles zero and TI rad. 

2 c 

Along each of these lines the Imaginary and real parts of Sq respec- 
tively are zero. We conclude that along the phase angle locus corr- 
esponding to zero rad., magnitude changes in L cause only phase changes 

in T. The reverse is true along the phase angle H rad. locus. 

2 

Since frequency response design work is generally carried out on 

the amplitude-phase plane, the mapping of the magnitude sensitivity 

function loci onto this space will now be considered. 

111-4 Mapping of the Magnitude Sensitivity Function Loci onto the 
Amplitude-Phase Plane . 

The magnitude sensitivity function loci can be trapped onto the 
amplitude-phase plane simply by inserting the variable change 

z - P 

into Equation (3-53). There results from this substitution, 



e ZZ + 2e Z cos 0 + 1 ■ ^ 1 j Z (3-61) 

Plotting this on the Z-0 plane gives the desired mapping. 

The ultimate purpose for the creation of this mapping is to use 
it as an overlay in design, just as is used in the Nichols Chart. By 
a simple manipulation it will be shown that the Nichols Chart itself 
is already calibrated for our use. In previous sections we concluded 
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„T o C 

S K “ S 0 



T+l" < 3 “ 62 > 

where, L(s) « G(s) H(s), Figure 3-5. It Is Intended to point out 
that this sensitivity function is obtainable directly from a convent- 
ional Nichols Chart, given the open-loop frequency locus 

L(ji» - G(jto) H(j») 

for the linear system in question. The procedures Is quite direct. 
Both numerator and denominator of Equation (3-62) may be divided by 
L(s) giving, 

r • r t fm\ * 

(3-63) 



S C - S T 
0 K 



LCs)~ 
1 + L(s) 



-1 



We seek the polar coordinate mapping of lines of constant magnitude 
sensitivity function given by 



1 


kl-kl 


- M 




or 


LCs)’ 1 


- M 


(3-64) 




1 + L(s)*’ 1 


let 


L(s) 1 • a 


i + jy 


(3-65) 


therefore 


* + jy 


H Vf 


(3-66) 




1 + x + jy 


* rl 



and from Equation (3-66) we get 




We recall the equation for the M circles of conventional design 
(ref. 9) and note that it is identical to Equation (3-67) . Thus we 
have the interesting result that the M-circles are to the open-loop 
frequency response curves as the magnitude sensitivity function loci 
are to the inverse open-loop frequency response. The Nichols Chart 
is just a plot of the M-circles (as well as the N-clrcles of constant 
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closed-loop phase) on the amplitude-phase plane, therefore the mapping 
of the magnitude sensitivity function loci onto the maplitude-phase 
plane is the Nichols Chart, provided that we read intersections as 
derived from inverse open-loop frequency response curves. 

The combination of the logarithmic amplitude scale and linear 
phase scale on the amplitude phase plane, with the symmetric Nichols 
chart is shown on Figure 3-9. Consider a typical open-loop locus and 
its inverse, on the extended Nichols chart as in Figure 3-9. The 

o | 

inverse function L ( jto) is obtained @n a point by point basis graph- 
ically and is symmetric with the direct function L(jfo) with respect 
to the point (0 DB, o deg). The Nichols chart is shown over a 340 
degree range and is seen to be symmetric (Ref. 9) about the axis 
phase-0 deg. 

One further simplification may be made as shown in Figure 3-10 
If we imagine that the phase calibration for L ‘(jk) locus on the 
Nichols Chart is opposite to that commonly used then the same M- 
contrours used in ordinary control system design, measure sensitivity 
function magnitude for use to the extent of being properly labelled. 

In this case L(jO) and L~*(jw) are symmetric with respect to zero DB 
line, with phase measured in opposite directions for each. For the 
particular example of Figure 3-10 we see that the maximum magnitude 
sensitivity function value is M ® 2.0. 

If it was desired, the magnitude sensitivity function for the 
closed loop system as a function of frequency could be read directly 
off the Nichols Chart as the intersections of the H-contours with 
the inverse open-loop frequency response curve. 
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The simplification Introduced by foregoing interpretation of 
the sensitivity function magnitude and phase loci could be very use- 
ful in the design of linear systems with regards to sensitivity 
specifications. 

Special Case : 

If the forward elements can be separated into two tandem parts, 

one of which displays sensitivity to the environment (denoted by an 

"environmental parameter" x), viz., 

G(s, x) - Gj^s) G 2 (s, x) (3-68) 

Then we may compute the sensitivity of the closed loop locus of the 

T 

overall system to the environmental parameter (denoted S^) according 
to 



or 



C dT/T dT/T 

l) " dL/L " dx/x 




o 



dx/x 

dL/L 

(3-69) 

(3-70) 



In the present example 

„o 



dL/L 

dx/x 



dG 2 (s,x)/G 2 (s,x) 

dx7x 



Therefore 



S T - 
x 



_C _ G 2 
S 0 S x 



(3-71) 

(3-72) 



where the first term of the product is evaluated graphically, and 

the second term is found by differentiation. It should be noted that 

if the parameter x appears in G 2 (s) as G 2 (s) ■ xG^(s), then 

G. ^ _ 

S « 1 and S - S G , which is in accordance with our relation 

x x 0 ’ 

C T 

S Q » , where k appears as in Figure Al-1 of Appendix 1. 
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Fig. 3-9: Typical L and its Inverse on Extended Nichols Chart. 
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III-5 - Philosophy of the Frequency Response Approach to Sensitivity 
Problem. 

Case 1 = When t . «* 0. 

ol 

In previous sections we discussed how sensitivity and frequency 
are related and how by utilising the Nichols Chart, the magnitude and 
phase of the sensitivity function and its peak magnitude can be obtained. 
Another approach which could be of high importance in most design cases 
is discussed as follows: 

Sensitivity specifications could be given in terns of the systems 

behavior on the jto axis. It is possible that the allowed variation 

in T(jto) may be exceedingly small over one frequency range and fairly 

large in other frequency ranges. The specification may obviously be 

in a wide variety of forms. For example the nominal T^, may be as 

shown in Figure 3-11, with bandwith iO, , amplitude peaking M , and the 

1 o 

specification may dictate that despite parameter variations, the 
magnitude peaking should never exceed 30 percent, the bandwidth should 
be in the range 0.8o^ to 1.2W^, and that from aero frequency to 
0.50^, | T| should never be less than 0.95. Finally, some thought 
should be given to the intermediate and the far off frequency ranges. 

For example, the behavior shown by dotted lines in Figure 3-11 is 
usually intolerable, as it results in high frequency ringing. A 
reasonable statement might be that for 1 . < to < , )T(jto)|< 3DB 

and for CJ > T(jW) must not exceed a predetermined value. 

The basic philosophy of the frequency response approach to 
sensitivity reduction is deduced from Equation 3-17, rewritten here 
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This equation applies whenever the 
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Fig. 3-13: To Obtain Permissible Location of -L (16J) 
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and therefore could be used in most control problems. Suppose that 
at some specific frequency s «* plant parameter variations are 

such that (k Q /kp at s ® may lie anywhere inside the region 
indicated in Figure 3-12. Suppose also that - ^(jA^) l 8 given by 
the complex number A + jB located at Q. Then Equation (3-73) states 

that T — » 

J_<L_ „ __9I_ 

T — * 

f QN 

The range of variation of the vector Q? fixes the range of variation 

of T^/T^ at s ** Suppose for example, that it is required that: 

1.20 > T /T c > 0.80 
o f 

at s » The problem is then to find the locus of -L q { which 

barely satisfies this inequality. A suitable procedure is as follows: 

I T I 

o I >0.80 with 
T f ' 

point N «= 1 + j as center, draw a circle of any radius R (Figure 
3-13). Next draw a circle Cj of radius 0.8R, using as center the 
point (the boundary of k Q /k^) which appears to be the closest to the 
first circle - this appears to be point J in Figure 3-13. The inter- 
section of the two circles determines a point Q. As a final check, 
use Q as center, and check whether a circle of radius 0.8R cuts k Q /kj 
locus only at J. If so Q is on the boundary of the locus of 
j T q /T gj > 0.80 at s ■ j^. The above construction is repeated, using 
different values for radius R, until the locus of Q is obtained. Such 
a locus for 1.20 >| T q /T^ | > 0.80 is shown in Figure 3-13. Clearly, 
mu6t be located inside the indicated region in Figure 3-13. 

On the other hand if there is the additional requirement that the 
phase of must not change more than 20°, then the permissible 

range of ”1^06^) reduced to the smaller region shown in Figure 3-13, 
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